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Summary. — The paper introduces a reader to a relatively young field of the 
physics of strongly interacting and disordered 2D electron system, in particular, to 
the phenomena of the metallic conduction and the apparent metal-insulator tran- 
sition in 2D. The paper briefly overviews the experimental data on the electron 
transport, magnetotransport, and spin-magnetization in 2D and on the electron- 
electron interaction effects. 
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1. Introduction 

Resistivity of materials around us ranges by more than 30 orders of magnitude, from 
10 20 to 10 -12 Ohmm. The variety of materials are classified as Metals (M) and Insulators 
(I), which are two distinct classes of matter only in the T — > limit. The principle 
difference between the two classes is in the character of the electronic wavefunctions, 
which are spatially localized in insulators and extended in metals. 

It is well-known, that subjected to even insignificant changes of pressure, dopant 
density, etc. many materials exhibit transformations between the two states P0|. The M- 
I transitions (MIT), unrelated to changes in the lattice structure and symmetry, are of a 
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special interest because they are considered to be continuous Quantum phase transitions, 
occurring at T = 0. The M-I transitions usually take place as density of electrons 
decreases; the latter necessarily leads to the increase of the effective strength of e — e 
interaction. A dimcnsionless ratio r s of the potential (Coulomb) energy to kinetic (Fermi) 
energy is often used to quantify interaction in 2D, r s oc n -1 / 2 with n being the electron 
density 2 . On the other hand, the decrease in the density is accompanied with the 
increase in the effective disorder, particularly, due to the weakening of the screening of 
potential fluctuations. Thus, both, interactions and disorder are important in the vicinity 
of the MIT; the separation and interpretation of their effects represents a hard task |3] . 

In the limit of strong interactions r s — > oo and zero disorder, the groundstate of the 
2D system is believed to be the Wigner crystal (WC) of electrons. Correspondingly, at 
zero temperature and zero disorder, for 2D system, there is at least one critical point on 
the interaction axes, quantum melting of the Wigner crystal 0J|2]- A weak disorder is 
expected to stabilize the WC El - As disorder increases further, WC looses long-range 
order and crosses over to a localized phase. Thus, in the limit of strong disorder and 
weak interactions, the ground state of the system is built of the single-particle localized 
states T . The experiments always take place at non-zero disorder and temperature. The 
subject of this paper is related to phenomena which occur due to the combined action of 
strong inter-electron interactions and disorder. 

The outline of this paper is as follows: Sections 1 and 2 introduce a reader to the 
field and briefly review the main results in the semiclassical and quantum frameworks, 
correspondingly. Section 3 describes major results on transport in 2D systems, which are 
discussed in terms of the apparent metal-insulator transition in 2D. Section 4 describes 
quantitative studies of the effects of electron-electron interaction; the experimental results 
are compared with available theories in various regimes. 

1.1. Mott semiclassical picture of the MIT. 
At finite temperature, or for finite system size L, the metal and insulator states are not 
well-defined; one can be mislead by attempting to distinguish these two states according 
to their conduction. In fact, the two states are inherently related with two distinct types 
of electron wavefunctions, localized and extended ones. The latter are classified according 
to the localization length £: if £ > (ra — r\) the waves are extended (metal) and \& cx 
cos(fcr); if £ < (r 2 — r\) the waves are localized (insulator) and \& cx exp(— r/£) cos(fcr). 
Note, the electron spacing, (r 2 — r{) ~ Xf — 27r/kp — 1120 (l0 11 .g l ,/2n) 1 ^ 2 A (the valley 
degeneracy g v = 2 corresponds to (lOO)Si-crystal plane). 

For the extended states, when Xf <C I < L (with I being the mean free path) the 
electron motion results in a classical Drude conduction a — ne 2 T p /m — ne 2 l/vF'm. As 
disorder increases, I decreases. However, / can not be less than Xf (the so called Ioffc- 
Regel criterium (I960)): Z m ; n = Xf = 2n/kF- Thus, we obtain the minimum metallic 
conductivity P CJ" min (Z m i n ) = 2ire 2 /h in 2D. We used here kp D = V2wn. Note, that 
h/e 2 — 25812.7 Ohm. In this semiclassical picture [Q, the M-I transition is caused by 
the disappearance of the localized (or, vice versa, extended) states. 
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1.2. Basics of the 2D semiconducting devices 
The 2D electron and hole systems discussed in this lecture consist of the bulk carriers 
confined in a two-dimensional potential well; the latter is formed at the interface between 
two semiconductors or semiconductor/insulator interface. Figure 2 schematically shows 
crossview of the silicon metal-oxide-semiconductor (MOS) stucture The 2D layer 
of electrons is confined at the interface between Si and SiC>2, grown at the top of the 
Si-substrate. On the top of the insulating layer, a thin metallic film (gate) is deposited. 
Throughout this lecture, the electrons and holes are the same quasiparticles and for 
simplicity will be often called "electrons" . 

When the positive voltage V g is applied to the gate (relative to the source, or drain 
contact, or relative to the Si-bulk), the conduction E c and valence E v band edges bend 
down as Figure [21 shows. For a sufficiently high gate voltage, the bottom of the con- 
duction band decreases below the bulk Fermi level and the resulting triangular potential 
well starts populating with electrons J2J. At room temperatures, when the bulk Si is 
conducting, redistribution of the electrons between the bulk and the 2D layer occurs in 
accord with the Poisson equation \. When the gate voltage is applied or changed at 
low temperatures, the bulk conduction is frozen-out, the equilibrium with bulk is not 
achieved, and the required electrons come into (or out) the 2D layer from the potential 
contacts. The contacts are lithographically defined dopant diffusion areas with high con- 
centration of electrons supplied by the dopants. For our purpose, it is essential only that 
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Fig. 1. - Schematic structure of a silicon field 
effect transistor. The 2D electron gas is con- 
fined at the Si side of the Si/SiC>2 interface. 
Drain and source are diffusion areas used as 
ohmic contacts to the 2D layer. Metallic gate 
is deposited on the top of the Si02 layer 
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Fig. 2. - Schematic energy band diagram of 
the Si-MOS structure. E c {Ev) is the bot- 
tom (top) of the conduction (valence) band, 
i?A-bulk acceptors energy level, Ef - Fermi 
level, Z-axes is directed to the bulk, perpen- 
dicular to the surface. 



the layer of electrons confined at the interface is two-dimensional, and that its charge Q 
is controlled by voltage V g applied to the gate: Q = ne = C(V g —Vt), with Vt, the so 
called threshold voltage. Thus, the Si-MOS system may be viewed as a plane capacitor, 
whose one plate is formed by the 2D layer of electrons whereas the gate serves as another 
plate. 
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2. Quantum transport at zero field 

2.1. Various types of transport: Delocalized states: diffusive and ballistic transport. 
There are two regimes of the metallic-like transport: diffusive for I < l v < L, and ballistic 
regime for I > L. The phase breaking length l v = ^JDt v (with t v , the phase breaking 
time, D - diffusion coefficient) is related with large changes of electron energy and is 
similar (not equivalent though) to inelastic length. The relevant scattering processes for 
t v at high temperatures are (i) electron-phonon scattering, and for low temperatures - 
(ii) electron-impurity and (iii) electron-electron scattering. In general, t v oc T~ p with 
p = 1 — 2, depending on the dominating scattering process. Since electron-electron 
collisions conserve the total momentum, the transport time l p (momentum relaxation) is 
determined by the two former processes and is insensitive to electron-electron scattering. 

2.2. Localized states . 

If an energy barrier A separates the energy of the electron states from conduction (or 
valence) band, and the wavefunctions of localized states don't overlap, the conduction oc- 
curs via temperature activated transitions whose probability is a oc exp(— A/T). When 
the probability of the temperature activation is too low (e.g., for high A, or for low 
temperatures), transport occurs via elastic tunneling between the localized states. This 
is the so-called hopping conduction regime, where the characteristic temperature depen- 
dence of the conductivity is a oc exp((— Tq/T) p ), where p = 1/2 or 1/3, depending on a 
specific model 7 ■ Usually, hopping conduction dominates at low temperatures. 

2.3. Electron's Phase Coherence and Transport. 

In the quantum-mechanical picture, the electron waves propagate and interfere. The 
interference gives rise to quantum corrections to the conductivity (SJ E| . In 2D system of 
non-interacting electrons (Fermi-gas), as T — > 0: 

2 2 

(1) o = (7jj — — ■ Wjv/t) ~<t d + ^— ln(Tr), 

7r/i 7rn. 

where ap is the semiclassical Drude conductivity. The single-particle interference thus 
may be viewed as a quantum "backscattering" . Weak logarithmic dependence is a typical 
attribute of the 2D system at high electron density (low r s ), and may be observed at 
high conductivity a 3> e 2 /h in the low-temperature diffusive regime Tr <C 1, r <S t v (see 

%B). 

2.4. Suppression of the weak localization in H± fields 

The quantum interference corrections originate from small loops of self-intersections on 
the trajectories, which quantum particles propagate in both directions, due to the time 
reversal symmetry 8 9. The amplitudes of the wavefunction for a particle passing the 
loop clockwise and counterclockwise acquires additional factors: A\ — > Ai exp (±z7r$/$o) 
with $o = h/e being the flux quantum and $ the magnetic flux through the loop. The 
magnetic field applied perpendicular to the 2D plane of motion, destroys the interference 
by reducing the probability for a particle to return. The interference breaks down when 
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Fig. 3. - u(T) dependences for high densities (from top to bottom, in units of 10 11 cm -2 ) 
32.1, 42.94, 48.4, 53.8, 64.7, 75.6, 86.5. 



the phase difference Aip becomes ~ 1: 

(2) A,>~^~1 H~^- 

The "negative magnetoresistance" effect is a tool to determine t v |9l 111). 

2.5. Single-particle scaling theory of localization (E. Abrahams et al., 1979) [T^] 
The one-parameter scaling theory considers the dependence of the conductivity on the 
system size. When the size L is changed, the effective disorder is changed; it is assumed 
that the only measure of this is the conductance G. The above assumption is equivalent 
to introducing a function (3 such that 

LdG d\nG 

[> 1 GdL dlnL' 

where (3 is a universal function of G solely. 

One can consider how the scaling hypothesis fits various obvious limits: 

(1) In the Ohm's law region (G 3> e 2 //i), (3 = d — 2 and the conductance does not depend 
on L for the two-dimensional case, d = 2. 

(2) If the states are localized, G(L) cx exp(-£/0- Therefore, (3{G) = \n(G/G ). 

The two distinct limits enable to plot the asymptotes of j3 in fig. 21a. One may expect 
that between these two limits, should vary smoothly, as shown schematically in fig.^la. 

Note: It is not necessary to vary L in order to move along the scaling trajectories. 
For finite temperatures, l v causes inelastic cut-off, so that the temperature tunes the 
effective sample length. 
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Fig. 4. - a) Schematic behavior of the scaling function for non-interacting electrons in different 
dimensions d. Bold dot depicts the critical point of the M-I transition for d — 3. b)Spin-orbit 
scattering case for d — 2. 

In d = 3, as fig. 01a shows, the scaling function j3 < for very small G (strong 
disorder). As conductance increases, passes through zero at a certain critical value G c . 
This is a repulsive point or an "unstable critical point" in the renorm-group terminology. 
For low initial G < G c , increasing the size of the system will cause G to decrease and 
the behavior of the system approaches that for the localized states. In contrast, when 
the initial conductance is high, G > G Cl increasing the size of the system will leads to 
the Ohmic regime (metallic conduction). Such behavior and the corresponding metal- 
insulator transition are consistent with the Mott' semiclassical picture described above. 

In lower dimensions, d < 3, (3 is always negative and all states are localized. As 
system size L increases, the system moves to the strongly localized regime. The scaling 
function describes a crossover from weak localization (localization length £ > L) to strong 
localization (£ < L) regime. 

It follows from the one-parameter scaling theory that "there is no true metallic con- 
duction in d < 3" [T2|. F° r d = 1, this conclusion is obvious since the negative scaling 
function has substantially large amplitude \(3\ > 1. In accordance with eq. (4), this causes 
the conductance to decrease quickly as L increases, approaching that for the strongly lo- 
calized case. However, the d = 2 case is more tricky: for high G >> 1, the scaling 
function is so close to zero, that minor correction may potentially change the sign of /3, 
therefore the situation requires a special attention. 

a) Spin-orbit scattering case (Hikami, Larkin, Nagaoka (1980)). The electron inter- 
ference is destroyed by spin-orbit random scattering and the weak localization "back- 
scattering" correction vanishes ^J|SJ|5]. More over, the spin-orbit scattering induces 
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a positive correction ("forward scattering") so that the 0- function becomes oc +1/2G 
for high G, as schematically shown in fig. 0]b. Note, that this results is valid only for a 
hypothetical 2D system of non-interacting particles. 

b) Interaction quantum corrections in the diffusive regime. Altshuler-Aronov correc- 
tions take into account e — e interactions in the Tr <C 1 limit [5| ^] : 



(4) 



5a 



2ir 2 h 



1 + 3 1- 



ln(l 



F" 



\n{h/k B TT), 



where Fq is the Fermi-liquid interaction parameter, dependent on r s (see further for 
more details). In the Fq —> limit, the total quantum correction <5o"i oc + Sa so + 8a ee 
is usually negative (backscattering). However, when Fq < and \Fq\ is large, the total 
quantum correction may lead to delocalization ( "forward scattering" ) . This case will be 
discussed in Sec. 0] in more details. 

c) Interaction quantum corrections to the transport in the ballistic regime Tt 3> 1 
have been calculated in Refs. ^3^)]; they are discussed in Sec. [1] 



3. An apparent MIT in 2D 

A great body of experimental data accumulated since 1979 support the conclusions 
of the scaling theory An example is shown in fig. 0a. The uppermost curves in 

figs. [5] show that the resistivity for a very disordered (low-mobility) Si-MOS sample at 
low densities exhibits an exponentially strong increase with cooling. This is consistent 
with motion to the left hand side along the logarithmic part of the /3-function in fig.0]a. 
For higher densities (lower curves in fig-JSJ), the exponential increase of the resistivity is 
replaced by a weak temperature dependence. The latter one at very high carrier densities 
in Si-MOS or in GaAs samples is logarithmic (see, e.g. fig. |3J and consistent with the 
theory {dp/dT > 0). 

As sample mobility increases, and, correspondingly, the density of the crossover (from 
the exponential to logarithmic T-dependence) decreases 1 the temperature derivative 
dp/dT obtains "wrong" positive sign (compare figs. [Ua, b, and fig. EJ. This effect was 
noticed more than two decades ago |2(JII21ll2"2"] and discussed in terms of the temperature- 
dependent screening 

|23 HUES]. As will be shown below, the sign of dp/dT is determined 
by the interaction parameter Fq, which in its turn, depends primarily on the electron 
density. 

The strong indication that the prediction of insulating states in 2D may not be uni- 
versally valid was obtained in 90s in studies of the quantum Hall liquid to insulator 
transitions |26l 1271 125 ] , Extended states, which in high magnetic field B are located 
at centers of Landau bands, were found experimentally to merge and remain in a finite 
energy range as B — > 0. This behavior is not expected within the scaling theory, which 
predicts that the extended states "float up" in energy in that limit 123 EH] ■ A serious 
challenge to the scenario of the insulating ground state in 2D arose in 1994 and the 
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Fig. 5. - Resistivity vs temperature for disordered Si-MOS samples: (a) with peak mobility 
fj, = 1500cm 2 /Vs, (b) fi = 5100cm 2 /Vs. The densities (in units of 10 11 cm -2 ) on the left panel 
span from 3.85 to 37.0 (from top to bottom), on the right panel - from 1.95 to 5.44. 

years following, when metallic behavior of electrical transport in zero magnetic field was 
observed in high mobility Si MOSFET structures [321 021 (see fig.®. 

Subsequently, qualitatively similar metallic states were found in other high mobil- 
ity 2D systems (for an extensive bibliography, see Ref. H2D- These include p— and 
n-Si/SiGe- [33 HH E7j, p-GaAs- EH1 EHl SOI ED E2, n-GaAs- 03 E], n-AlAs- 
heterostructures 05] i inverted Si-on-insulator structures an d back-gate biased Si- 
MOS structures |S7|. The most important findings of these studies can be summarized 
as follows: 

Features in zero magnetic field: 

• If the density of 2D carriers n is larger than a sample-dependent value n c , the p(T) 
dependence is metallic-like, i.e. dp/dT > (see fig.®. The changes in p can be 
almost an order of magnitude (in high mobility Si-MOSFETs |32|). 

• The strong metallic-like /d(T)-dependence is characteristic of a wide range of den- 
sities rather than of the critical regime; it sets in for relatively high temperatures 
T < Tp in the ballistic rather than diffusive regime ( see fig- Q • 

• When n < n c , both the sign (dp/dT < 0) and the functional (exponential) form of 
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7(K) 

Fig. 6. - Resistivity vs temperature for HIGH 
mobility Si-MOS sample 51 in the critical 
regime n ~ n c . Electron density decreases 
from 1.326 in steps of 0.0218 (from bottom to 
top) and in unites of 10 11 cm~ 2 . 




Fig. 7. - Resistivity vs T/Ef in a wide range 
of densities. The diffusive/ballistic border for 
Tt ^] is calculated based on the data [48] 
with renormalized m*. 



the p(T) dependence are characteristic of the insulating behavior (see figs. I5I6|) . 

• At much higher densities (n > 30 x n c ), the derivative dp(T)/dT is metallic (posi- 
tive) at higher temperatures T ~ Tp, but p exhibits an slow (logarithmic) "insulat- 
ing" upturn at lower T, consistent with eq. (2) ^J^HSU] - see fig.Ol Qualitatively, 
conduction at high density is consistent with the weak localization picture. 

• The log(i?(T)) data show a reflection symmetry with respect to n c |.i2l I33| (see 
fig. EJ), with exclusion of the low-temperature data, which in the T — ► limit 
diverge for n < n c , but do not fall to zero for n > n c . This signifies that at 
intermediate temperatures, the p(T)-dependence is empirically represented by a 
logarithmic function p(T) oc exp [±(Tq/T) p ] with p ~ 1 and the sign equal to that 
of (n c — n) |51j ; such dependence obviously fits the one-parameter scaling function 
in the vicinity of a critical point shown in fig. 0]b. 

• For high-mobility samples, the resistances of the metallic and insulating phases, 
measured as a functions of T and n, exhibit scaling properties with respect to 
T/Tq [62i 1331 loT] demonstrated in fig. [S] The energy scale To goes to zero on 
both sides of the transition; such behaviors are reminiscent of a critical point of a 
continuous quantum phase transition |52| . This interpretation, however, seems to 
work only for intermediate temperatures and fails in the T — > limit. For a more 
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detailed review and phenomenological 
critical regime, see ref. 
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Fig. 8. - Resistivity for high-/i sample scaled 
vs T /T [5T] 



discussion of the experimental data in the 
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Fig. 9. - R(H«) dependences for different den- 
sities |59|. 



Features in the in-plane field 

• In Si-MOSFET and p-GaAs structures, p increases strongly with H\\ [5JJ |55] EHl 
EZIEHHSniEni and shows saturation at a field, which is approximately equal to that 
for complete spin polarization (HJ; this is most sharply pronounced in Si-MOS 
structures |551 1551 151155] . 

• As B\\ couples to spins rather than to orbital motion, it follows that the spin-related 
mechanism of interactions (i.e. exchange) is of major importance. More specifi- 
cally, when ill I field switches off the spin-degrees of freedom, the R(T) dependence 
flattens and MIT disappears (S21ES]- 

As n c is quite small (~ 10 10 — 10 11 cm" 2 ), the corresponding r s values are large, 
an order of 10; therefore, it is reasonable to assume the e — e interaction to be one of 
the major driving forces in the above listed phenomena. In the following section, the 
e — e interaction effects will be discussed in more detail. We shall show, in particular, 
how the listed features find a natural explanation within the frameworks of the quantum 
interaction corrections at intermediate temperatures. 

4. Quantitative studies of the electron-electron interactions 

4.1. Fermi- liquid rcnormalization of electron parameters in 2D systems 

The Landau Fermi- liquid (FL) theory |64l Ifi5| introduces a number of interaction param- 
eters to describe the interacting system in terms of the renormalized <?*-factor Lande, 
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Fig. 10. - Renormalized spin susceptibility (a), effective mass (b), and gr*-factor (c) measured 
with Si-MOS samples |ESII77| . 



effective mass m* etc.: 

(5) ^o=^-l, F > = 2 (£-1 

9 \m b 

Recently, these parameters have been determined experimentally in a wide range of 
densities, in measurements of the quantum oscillations of conductivity (Shubnikov-de 
Haas (SdH) effect) in tilted or crossed weak magnetic fields |^| 1661 1671 16"5| . In these 
experiments, the renormalized spin susceptibility \* is obtained with no assumptions: 

(6) X = 9 9bm fiB-r, 



Xb 2m b (F a + 1)' 

where g b , m b , and Xb are the corresponding band values. Figure 1101 shows the r s - 
dependence of the measured renormalized x* and m*-values |57l I66j : the renormalized 
(^-values are obtained from these two sets of data as x* / m * ■ Other experimental ap- 
proaches consist in scaling the a(H/ H )-d&t& for different densities HHH E3 E] > or m 
fitting the <r(T) and (t(Hh) dependences to the calculated quantum corrections \R>\, us- 
ing Fq (and in some cases, also m*) as fitting parameter |72l 1731 ITU I75| . 

The resulting Fq values obtained for 2D electrons in GaAs jHHl and in Si-MOS samples 
[7TlCnilMHZlC3IZilZ21 qualitatively agree with each other [77j but differ substantially 
from the data for 2D holes in GaAs (73], and, especially, from the data taken for low 
carrier density 021; the reason for this discrepancy is not clear yet |77| . 
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Figure lTOl shows a strong increase in m*, and g* with r s , which should significantly 
affect transport at low densities. The growth in \* might reflect a tendency to either 
ferromagnetic or antiferromagnetic transition; besides, within the FL-description eq. 
the gf-factor should diverge as F§ — ► — 1. Therefore, a separate, interesting question is 
whether or not m* and g* diverge as r s increases, and whether the divergency (if 
any) occurs at the MIT. Numerical calculations can not give an unambiguous answer 
[751 1701 150) for the single- valley system, whereas the 2D valley system is expected to be 
more stable and to remain unpolarized up to the WC phase. 

As a result of the renormalization, the energy scale Ep diminishes. However, in the 
vicinity of the transition (i.e. at the sample-dependent density n = n c ), m* does not 
diverge [51] [77]; this may be easily seen, e.g., from non- vanishing amplitude of the SdH 
oscillations measured at T = 0.03 K for r s values up to 9.5 [5T]. Thus, the Mi-transition 
is not caused by simply vanishing the Ep value and the quasiparticles. In Refs. (7TJ 100] 
the in-plane magnetoresistance p(-H)i) data have been discussed in terms of a developing 
ferromagnetic instability at or very close to n c . However, direct measurements of the 
(a) magnetization of 2D electrons (52], (b) period, and (c) sign of the SdH oscillations 
|811 177] do not confirm such possibility: the spin susceptibility remains finite at n > 
7.7 x 10 cm , the density which is lower than the critical n c - value for many samples. 

The strong enhancement of the weak field spin susceptibility was predicted also by the 
renorm-group (RG) theory [1411831 151] : however the RG-theory is based on consideration 
of the diffusive modes of electron-electron interaction in the regime Tt <C 1, whereas 
the renormalized parameters (Fig. 10) are measured in the ballistic interaction regime 
Tt > 1 (for a more detailed discussion, see Section 4.3). 

4.2. Implementation of the measured FL parameters to the metallic-like transport. 

A considerable progress has been also achieved in the theory of quantum transport: 
in Refs. ^3 E] the interaction corrections to the conductivity have been calculated 
beyond the diffusive regime, in terms of the FL interaction parameters. In this section, 
we make a comparison of the data with the theory and conclude that the "metallic" 
drop of conductivity with cooling in the regime a 3> e 2 /h can be accounted for by the 
interaction effects in electron "liquid" at temperatures which correspond to the ballistic 
regime Tt > 1. This observation suggests that the anomalous "metallic" conduction in 
2D, at least for densities not too close to the critical density, is the finite-temperature 
phenomenon rather than the signature of a new quantum ground state. 

The theory [Rj\ considers backscattering of electrons at the scattering centers and at 
the Friedel oscillations of the density of surrounding electrons. The interference between 
these scattering processes gives rise to the quantum corrections, which are calculated to 
higher orders in the interactions and leading order to the temperature. The interference 
gives rise to quantum corrections to the Drude conductivity (in units of e 2 /irfi): 



(7) 



cr(T) -a D = Sa c {T) + 15Sa T (T). 



Metal-insulator transition in 2D 



13 



Here 



(8) 



Sa c = 
Sa T = A(Fg)x 



i-lt(x,F a : 



--In 1^1 



2tt 



T 



2tt 



and 



E F 



T 



are the interaction contributions in the charge" (exchange term and singlet part of the 
correlation terms) and triplet channels, respectively; x — TTk B /h, A(Fq) — Fq/(1 + 
F a ), and B(Ffi) = 1 -ln(l + F a )/F a . The prefactor 15 to Sa T reflects enhanced number 
of triplet components: two valleys and two spins produce 4 pseudospin components; 
interaction of two particles involves 4 x 4 = 16 channels, one of them being singlet. 

Diffusive regime, Tr <C 1. For large \Fq \ and Fq < 0, 8cr(T) becomes positive. E.g., 

w 



for F a = -0.3 (i.e. g 



(9) 



15 x 0.2 



l)lnT = InT. 

' 2tt 



Thus, due to valley degeneracy, the quantum correction is constructive and produces a 
"forward scattering" . 



(10) 



Ballistic regime, Tt > 1. The quantum correction Sa(T) is quasilinear in T: 

pa 



5a(T) 



1 + 15, 



1 



(Tr) 



Comparison of the theory with experiment at high G ^> e 2 /h (see Fig. 

The terms in eq. JJJ are functions of Tt/Ti and F§. The momentum relaxation time 
r may be found from the Drude resistivity pr> = o - ^ 1 using the renormalized effective 
mass m* determined in Ref. |66| . Thus, the comparison shown in fig. II II has no fitting 
parameters, beyond an. The determination of the latter value is transparent |73 |. it 
is obtained by a quasilinear extrapolation of the experimental a(T) data to T = 0, by 
setting the logarithmic (diffusive) terms to zero. Almost entire temperature range in 
fig. ^2 belongs to the ballistic regime. In view of the quantitative agreement of the 
experimental data with theory, we conclude that the initial slope of the a(T) dependence 
is well described by the quantum interaction corrections eqs. Q with F£ (r s ) and m*(r s ) 
values determined independently (see fig. 1 1 0|) . Several experimental teams came to a 
similar conclusion for p-GaAs/AlGaAs [23 E2 and Si MOSFETs [7117511551173] . 

We conclude this section by listing the achieved results: 

• Not too close to the critical regime, a ^> e 2 /h, the "metallic" strong T-dependence 
results from the quantum interactions at intermediate range of T"s. 

• As T decreases further, in the diffusive regime, a(T) is expected to vary cx InT 
with da /(IT < for sufficiently large g* (large \Fq |) and large number of valleys. 
Then, for T < Tq (determined by the intervalley h/r vv and spin-spin scattering 
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Fig. 11. - Conductivity for sample Si22 vs temperature [85) . The electron densities, from top 
to bottom are: n = 5.7, 6.3, 6.9, 7.5, 8.1, 8.7, 9.3, 10.5, 11.7, 12.9, 14.1, 16.5, 18.9, 21.3, 23.7, 
28.5, 35.7 (in units of 10 11 cm -2 ). Dots represent the data, solid lines - the theoretical curves 
with F a and m* from Ref. gSl- 

rate h/r ss ), the enhancement factor (15) should vanish and the localizing terms 
are expected to overwhelm. 

• If the true metallic phase in 2D exists, it must of a non-FL type |22 E3 EH- 
Several non-FL theoretical models have been recently suggested, such as phase 
separation model (Wigner solid inclusions in the FL) an electron-paired state 
with broken U(l) symmetry |90| . magnetic phase transition etc. However, so far 
no clear signatures of the non-FL behavior have been detected experimentally. 

4.3. Transport in the critical regime, a ~ e 2 /h and n ~ n c 

In view of the encouraging comparison with the theory of interaction corrections in 
the high density regime, it seems attractive to extend the comparison to the critical 
regime, where a ~ e 2 /h (see fig. 01. Quantum corrections theory is inapplicable for large 
conductivity changes Safer ~ 1 and in the vicinity of the critical point, where a ~ e /h. 
The method now commonly in use is a generalization of the nonlinear cr-model theory, 
which has been developed by Finkel'stein ^3] and by Castellani, Di Castro et al. |83U91| . 
The renormalization-group method is an extension of the same diagrams of the interaction 
theory to the case of the strong disorder and interaction 92 . The RG equations describe 
rcnormalization of both, disorder (i.e. p) and interaction parameters as the length scale 
varies with T [HI I9T1 155] : 
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Fig. 12. - a) Comparison of the p(T) data normalized to its maximum value (1.35ft/e 2 ) with 
solution of the RG equations (solid line) eq. Hill : b) Renormalized interaction parameter 72 
|94|. calculated from the 2nd Eq.(ll) with p(T) data shown on the left panel. The „0"-point at 
the ln(TY) scale corresponds to T = 1.8K. Sample Sil5, electron density n — 0.88 x 10 11 cm -2 . 



(11) 



d P P 2 \^ 1 1 q ( 1 + ^ 2 wi , ^ 1 
— = — 1 + 1-3 ln(l + 72) - 1 

»!/ t [ \ 72 

^72 _ pQ- +I2) 2 
dy 2n 



where y = — ln(Tr), 72 is the Fermi- liquid amplitude for large angle scattering (72 — > 
—Fq/2 in the weak-coupling limit), and p is in units of /1/e 2 . The 1st term in the square 
brackets is the weak localization contribution, the 2nd term is the e — e interaction 
correction in the singlet channel ,8, 9 , and the last term is the contribution from 3 triplet 
modes. Due to the difference in symmetry of the singlet and triplet wavefunctions, the 
singlet and triplet terms have different sign, and favors localization and delocalization, 
correspondingly. For two- valley system, factor 3 should be again replaced by 15, and the 
weak-localization term also becomes twice larger [M] . 

The above equations have a universal solution, which has been compared with experi- 
mental p{T) data in Ref. |23J. We present in fig.^|a similar comparison, which confirms 
a certain similarity between the p(T)-data and the theory. There is a limitation though, 
the eqs. (|11[1 are perturbative, i.e. are derived in the lowest order in p, therefore, the 
comparison with experimental data for p(e 2 /h) ~ 1 can be only qualitative. 

The application of the RG perturbative equations to the experimental data is based 
on the following assumptions: 
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(i) the critical regime n « n c belongs to the diffusive domain of e-e-interactions TV <C 1; 

(ii) the p{T) maximum signifies a turning point from localization to derealization be- 
havior due to the renormalization of the interaction parameter (primarily, 72) with T; 

(iii) the decrease in p with cooling takes place mainly for the account of the strong growth 
in the interaction parameter 72 (and hence, F§ ). 

The first assumption does not agree well with the empirical diagram in fig. El therefore, 
its application would require to re-interpret the definition of the diffusive/ballistic border, 
i.e. to find a different interpretation for the r(n) or Fq(h) data in the critical regime. 
In Ref. |94) . the assumption (ii) was presumed to be fulfilled in the most clean samples; 
however, some high mobility Si-MOS samples do not exhibit a pronounced p(T) maximum 
in the critical regime 123031 (see, e.g., fig.HjJ). Furthermore, the p(T) data in the critical 
regime, as a rule, show a strong sample-dependent and non- universal behavior |511 I9b| . 
It is possible though that the relevant disorder is related with short-range fluctuating 
potential at the interface which may favor inter-valley scattering and, thus reduce a 
number of the triplet terms in eq. (|llf> but is not significant for the mobility at low T. 

Regarding assumption (iii), figure IT2T b shows that the interaction parameter 72, deter- 
mined from the experimental data, starts growing at rather high temperatures, ~ 1.8 K 
for n Ri n c ; the latter behavior of 72 is necessary in the RG-theory to turn down the 
resistivity flow. Thus, the renormalization of 72 should occur in a quite accessible range 
of temperatures (0.3 - 1.8)K. However, the anticipated strong T-dependence of the renor- 
malized spin susceptibility, eq. 10, 83 is not observed experimentally at these temper- 
atures in SdH-effect measurements: x* /Xb oc g*m* /2rrib changes only within 2% in the 
above temperature range |77) . 

An alternative theoretical approach (S5| suggests that the critical behavior in the 
vicinity of n c is dominated by the physics of the insulating state, namely by the melting 
of the disordered Wigner crystal electron glass [§7|, or an inhomogeneous state 
consisting of inclusions of WC into the 2D liquid 89 . Corresponding theories has not 
reached yet a fully predictive stage. We note that signatures of the glassy behavior were 
reported in Ref. [5*5] . 

4.4. Transport in the presence of the in-plane field 
Figure El shows a typical behavior of p(B\\). At low densities n ~ n c , as field increases, 
transport becomes temperature activated (hopping regime) |63| . When n approaches 
n c , by definition, n c (H\\) tends to n c (Hu = 0) = n c , therefore, the range of magnetic 
fields where the magnetoresistance (MR) may be studied in the regime of "metallic" 
non-activated conduction, shrinks to zero. 

Regime of high densities n 3> n c , G ^ e 2 /h. According to the theory of interaction 
corrections 1151. 



where K b » {E z / 2T) 2 f {Fg) in the low field limit E Z /2T < 1, and K b » (E z /T)f 2 {F$) 



(12) 
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in the limit Ez/2T 3> 1 limit. Correspondingly, as field increases, Sa(H) should increase 
initially oc H 2 /T, and than oc H . 

Experimental data, in general, show similar behavior |86l 1661 I73| with a transition 
from parabolic to linear dependence; however, in contrast to the temperature dependence 
of p, the agreement with theory ^H] is only qualitative I n high fields, g*pH\\ ^> T, 
the magnetoresistance deviates substantially from the theory, and the deviations arc 
sample dependent [EH- As density decreases and field increases, the deviations of the 
measured MR from the theory grow and reach a factor of two at n = 1.5n c [731 I63| . 
The measured a{H,T) scales with temperature somewhat different from eq. (|12|) . which 
predicts 8a oc — (H 2 /T) in low fields. 

Critical regime of densities n ss n c , and a ~ e 2 /h. As mentioned above, the magne- 
toresistance studies in the critical regime are restricted to the fields vanishing to zero 
as n — > n c . Nevertheless, down to densities n ~ l.ln c . the magnetotransport can be 
safely studied in the regime g*(iH^ <C T. The experimental data in low fields scales as 
5a oc — (H 2 /T p ), where p increases from 1.1 to 1.6 as density decreases from 5n c to 1.2n c . 
Comparing this empirical scaling law with that predicted by the RG theory |141 , 

e 2 

(13) a(H, T) - a(0, T) = -0.084— 72 ( 72 + 1) 

7T/1 

we conclude that the interaction parameter 72 (roughly, oc \F§\) decreases as temperature 
decreases. This result is not consistent with the main idea of the two-parameter scaling, 
where 72 is expected to increase and x* diverges oc T -4 / 3 [S3] as temperature decreases 
(i.e., the length scale increases). It is important also that the direct SdH measurements 
in low if 11 fields do not confirm strong T-dependence of Fq |77| . It might be that 
the critical regime where the RG-description is applicable is much narrower, than the 
range of densities (« 20%n c ) where the p(T) data exhibit a pronounced maximum. If 
this is the case, the p(H\\ , T) analysis should be restricted also to much lower magnetic 
fields. We note that right at the transition and for lower densities n < n c , the p{T) 
displays an approximately T-activated exponential dependence p{H,T) oc exp(A/T) 
with A oc (n — n c {H)) oc H [S3]; the resulting dependence exp(H/T) is not easy to 
distinguish from the (i//T) 2 -dependence. 

In refs. [53] [Ml EMI E3 , it was suggested that transport in the critical regime is driven 
not only by "universal" effects of interactions among the itinerant electrons, but also by 
interactions of itinerant electrons with localized ones, the latter issue (short-range aspect 
of interactions) is missing in most of the existing theories. 

Homogeneity There is an important issue, whether or not a spatial inhomogeneity 
develops in the 2D electron system as the carrier density decreases; such inhomogeneity 
was suggested to cause the percolation-type Mi-transition |1001 HOlj . Though this issue 
was considered in a number of theoretical papers, there is no clear experimental evi- 
dence for such mechanism to be the major driving force in the experimentally observed 
phenomena in high mobility samples. 
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5. Summary 

I will present below a brief pedagogical summary of the issues which has been learnt 
in the considered field, and those which still remain unanswered. 

• For high densities and high conductivities, in the ballistic regime Tt 3> 1, it is 
now widely accepted that the metallic T-dependence of the conductivity is a fi- 
nite temperature effect caused by e — e interactions. However, (i) in the diffusive 
regime Tt « 1, a thorough comparison with theory is missing; some data reveal 
a substantial disagreement [73| with theory, which may be attributed, e.g., to the 
intervalley scattering; (ii) a successful comparison with interaction theory has not 
been demonstrated yet for the same set of p(H, T)-data for both zero and non-zero 
in-plane field. 

• For low densities, in the critical regime n « n c , the observed one-parameter critical 
behavior of conduction a(T) seems to be characteristic of intermediate tempera- 
tures only. For much lower temperatures, in the diffusive regime, the potential 
critical behavior remains unexplored. 

• So far, all data find a reasonable explanation within the FL theory; no clear signa- 
ture of the non-FL behavior was detected. 

• Both, the growth of the spin susceptibility as density decreases and the non- 
monotonic T-dependence of resistivity in the critical regime n ss n c , are in a qual- 
itative agreement with the RG theory. However, some other experimental results 
are inconsistent with the theory; the discrepancy requires more detailed analysis of 
the experimental data and the RG scenario of the MIT in 2D. 

• The behavior of the renormalized spin susceptibility and effective mass for 2D holes 
at high r s values remains puzzling; the reason of the discrepancy between the 2D 
hole and electron systems is not clear. 

• By now, there is no complete theory which would include short-range and high- 
energy aspects of e — e interactions, inter-valley scattering, etc. 

• The interesting ideas of the potential spontaneous polarization transition in the 
spin or valley system require investigations at substantially higher r s values. 

• The detailed experimental verification of the possible formation of local moments, 
melting of the Wigner glass and other non-FL scenarios is currently lacking. 
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